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A POSTERIORI ERROR ESTIMATES FOR FINITE ELEMENT 
EXTERIOR CALCULUS: THE DE RHAM COMPLEX 

ALAN DEMLOW^ AND ANIL N. HIRANI^ 



>»^ Abstract. Finite element exterior calculus (FEEC) has been developed over 

vN the past decade as a framework for constructing and analyzing stable and 

f\ accurate numerical methods for partial differential equations by employing 

^ differential complexes. The recent work of Arnold, Falk and Winther [4] in- 

^\ eludes a well-developed theory of finite element methods for Hodge Laplace 

^^ problems, including a priori error estimates. In this work we focus on develop- 

00 ing a posteriori error estimates in which the computational error is bounded by 
1— H some computable functional of the discrete solution and problem data. More 

precisely, we prove a posteriori error estimates of residual type for Arnold- Falk- 

1 I Winther mixed finite element methods for Hodge-de Rham Laplace problems, 
■^r While a number of previous works consider a posteriori error estimation for 
^^^ Maxwell's equations and mixed formulations of the scalar Laplacian, the ap- 
^; proach we take is distinguished by unified treatment of the various Hodge 

r^ Laplace problems arising in the de Rham complex, consistent use of the lan- 

guage and analytical framework of differential forms, and the development of 
a posteriori error estimates for harmonic forms and the effects of their approx- 
imation on the resulting numerical method for the Hodge Laplacian. 



> 

(^ 1. Introduction 



In this paper we study a posteriori error estimation for finite clement methods 
for the Hodge Laplacian for the de Rham complex generated by the Finite Element 
\IJ Exterior Calculus (FEEC) framework of Arnold, Falk, and Winther (abbreviated 

^^ AFW below). Finite element exterior calculus has been developed over the past 

^_H decade as a general framework for constructing and analyzing mixed finite element 

;;■ methods for approximately solving partial differential equations. In mixed meth- 

ods two or more variables are approximated simultaneously, for example, stresses 
and displacements in elasticity or pressures and velocities in fluid problems. The 
essential feature of FEEC is that differential complexes are systematically used in 
order to develop and analyze stable and efficient numerical methods. Historically 
speaking, some aspects of mixed finite element theory such as the so-called "com- 
muting diagram property" (cf. [10]) are related to differential complexes, and some 
early work by geometers such as Dodziuk [15] and computational electromagnet- 
ics researchers such as Bossavit and others [7] also contains ideas related to finite 
element exterior calculus. However, around 2000 researchers working especially in 
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2 A .DEMLOW AND A.N. HIRANI 

electromagnetics and elasticity [16, 2] independently began to realize that differen- 
tial complexes can be systematically exploited in the numerical analysis of PDEs. 
This work has culminated in the recent publication of the seminal work of Arnold, 
Falk, and Winther [4] containing a general framework for FEEC (cf. also [3]). 

Error analysis of numerical methods for PDE is generally divided into two cate- 
gories, a priori and a posteriori. To fix thoughts, let u solve — Au := / in a polygonal 
domain fl with Neumann boundary conditions ^ = on dil and J„ u = assumed 
in order to guarantee uniqueness. Also let u^ (z Sh he a finite element approxi- 
mation to u, where Sh C H^{n) is the continuous piecewise polynomials of fixed 
degree r with respect to a mesh Th- A classical a priori estimate is 

(1.1) \\u-Uh\\m<Ch-\\u\\r+l. 

Such estimates are useful for verifying optimality of methods with respect to poly- 
nomial degree and are commonly used to verify code correctness. However, they 
provide no information about the actual size of the computational error in any given 
practical problem and often assume unrealistic regularity of the unknown solution. 
A posteriori error estimates provide a complementary error analysis in which the 
error is bounded by a computable functional of u^ and /: 

(1-2) \\u-Uh\\<£{uh,.f). 

Such estimates provide no immediate information about asymptotic error decrease, 
but do ideally yield concrete and reliable information about the actual size of the 
error in computations. In addition, £{uii,f) and related quantities are typically 
used to derive adaptive finite element methods in which information from a given 
computation is used to selectively refine mesh elements in order to yield a more 
efficient approximation. We do not directly study adaptivity here. 

While there are many types of a posteriori error estimators [1, 5], we focus our 
attention on residual-type error estimators. Roughly speaking, residual estimators 
are designed to control u — u^ hy controlling the residual / + Au^ , which is not a 
function (since Vuh is only piecewise continuous) but is a functional lying in the 
dual space of H^{n)/M.. Given a triangle X € Th, let hx = diam(r). We define 
the elementwise a posteriori error indicator 

(1.3) 7j{K) - hK\\f + AuhU,^K) + h]['\\lVuhj\\L,(aKy 

The volumetric residual /i/cH/ + Aufi\\]^^tj^\ may roughly be seen as bounding the 
regular portion of the residual / + Au^ . jVu^J is defined as the jump in the normal 
component of \/uh across interior element boundaries and as Vm^ • n on element 
faces e C dfl. Since natural boundary conditions are satisfied only approximately in 
the finite element method, the latter quantity is not generally 0. The corresponding 
term in (1.3) may be thought of as measuring the singular portion of the distribution 
/ -|- Au/j. A standard result is that under appropriate assumptions on Th, 

(1-4) \\u - Uh\\mm/M <c{J2 viKff'^. 

Ken 

That is, £{ufi, /) ~ C^i^KeT vi^)"^)^ is ^ reliable error estimator for the energy 
error \\u — Uft||_f/i(a)/R- An error estimator £ is said to be efficient if £{uh,f) < 
C\\u~Uh\\, perhaps up to higher-order terms. Given K € Th, let ojk be the "patch" 
of elements touching K. We also define osc{K) = hxWf — -P/IUaCif)' where Pf is 
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the L2{K) projection onto the polynomials having degree one less than the finite 
element space. Then 

(1.5) n{Kf<ci\\u-ua]i.(^,)+ Y. osc(X)M. 

In our development below we recover (1.4) and (1.5) and also develop similar results 
for other Hodge Laplace problems such as the vector Laplacian. 

We pause to remark that residual estimators are usually relatively rough estima- 
tors in the sense that the ratio £{uh, f)/\\u — Uh\\ is often not close to 1 as would 
be ideal, and there are usually unknown constants in the upper bounds. However, 
they have a structure closely related to the PDE being studied, generally provide 
unconditionally reliable error estimates up to constants, and can be used as build- 
ing blocks in the construction and analysis of sharper error estimators. Thus they 
are studied widely and often used in practice. 

In this work we prove a posteriori error estimates for mixed finite element meth- 
ods for the Hodge Laplacian for the de Rham complex. Let HhP — >■ iJA^ ^ • • • ^ 
ff/^n-i _^ j^^ |-|g -j-j^g n-dimensional de Rham complex. Here A*^ consists of fc-forms 
and HK*' consists of L2-integrable fc-forms a; with L2 integrable exterior derivative 
duj. For n = 3, the de Rham complex is 7J^ — >■ iJ(curl) — >■ iJ(div) — >■ ^2- For 
Q <k <n, the Hodge Laplacian problem is given by 5du + dSu = f, where 6 is the 
adjoint (codifferential) of the exterior derivative d. When n = 3, the 0-Hodge Lapla- 
cian is the standard scalar Laplacian, and the AFW mixed formulation reduces to 
the standard weak formulation of the Laplacian with natural Neumann boundary 
conditions. The 1- and 2-Hodge Laplacians are instances of the vector Laplacian 
curl curl —Vdiv with different boundary conditions, and the corresponding FEEC 
approximations are mixed approximations to these problems. The 3-Hodge Lapla- 
cian is again the scalar Laplacian, but the AFW mixed finite element method now 
coincides with a standard mixed finite element method such as the Raviart-Thomas 
formulation, and Dirichlet boundary conditions are natural. We also consider es- 
sential boundary conditions below. 

Next we briefly outline the scope of our results and compare them with previous 
work. First, in the context of mixed methods for the scalar Laplacian and espe- 
cially FEM for Maxwell's equations two technical tools have proved essential for 
establishing a posteriori error estimates. These are regular decompositions [17, 22] 
and locally bounded commuting quasi-interpolants [24]. Relying on recent analyti- 
cal literature and modifying existing results to meet our needs, we provide versions 
of these tools for differential forms in arbitrary space dimension. Next, our goal 
is to prove a posteriori estimates simultaneously for mixed approximations to all 
fc-Hodge Laplacians (0 < fc < n) in the de Rham complex. Focusing individually 
on the various Hodge Laplace operators, we are unaware of previous proofs of a 
posteriori estimates for the vector Laplacian, although a posteriori estimates for 
Maxwell's equations are well-represented in the literature; [6, 24, 27] among many 
others. The estimators that we develop for the standard mixed formulation for the 
well-studied case of the scalar Laplacian are also modestly different from those pre- 
viously appearing in the literature (cf. §6.4 below). In addition, our work extends 
beyond the two- and three-dimensional setting assumed in these previous works. 
Throughout the paper we also almost exclusively use the notation and language 



4 A .DEMLOW AND A.N. HIRANI 

of, and analytical results for, differential forms. The only exception is §6, where 
we use standard notation to write down our results for all four three-dimensional 
Hodge-Laplace operators. This use of differential forms enables us to systematically 
highlight properties of finite element approximations to Hodge Laplace problems 
in a unified fashion. A final unique feature of our development is our treatment of 
harmonic forms. In §2.4, we give an abstract framework for bounding a posteriori 
the gap between the spaces S)^ of continuous forms and i^^ of discrete harmonic 
forms (defined below) . This framework is an important part of our theory for the 
Hodge Laplacian and is also potentially of independent interest in situations where 
harmonic forms are a particular focus. Since our results include bounds for the 
error in approximating harmonic forms, our estimators also place no restrictions 
on domain topology. We are not aware of previous works where either errors in 
approximating harmonic forms or the effects of such errors on the approximation 
of related PDE are analyzed a posteriori. 

We next briefly describe an interesting feature of our results. The AFW mixed 
method for the fc-Hodge Laplace problem simultaneously approximates the solu- 
tion ?i, (T ~ Su, and the projection p of / onto the harmonic forms by a discrete 
triple (ahjUfnPh)- The natural starting point for error analysis is to bound the 
HK^~^ X HK^ X L2 norm of the error, since this is the variational norm naturally 
related to the "inf-sup" condition used to establish stability for the weak mixed 
formulation. Abstract a priori bounds for this quantity are given in Theorem 3.9 of 
[4] (cf. (2.8) below), and we carry out a posteriori error analysis only in this natural 
mixed variational norm. Aside from its natural connection with the mixed varia- 
tional structure, this norm yields control of the error in approximating the Hodge 
decomposition of the data / when 1 < fc < n — 1, which may be advantageous. 

As in the a priori error analysis, the natural variational norm has some dis- 
advantages. Recall that residual estimators for the scalar Laplacian bound the 
residual / -I- Am^ in a negative-order Sobolev norm. For approximations of the 
vector Laplacian, establishing efficient and reliable a posteriori estimators in the 
natural norm requires that different portions of the Hodge decomposition of the 
residual / — dah — Ph — 6duh be measured in different norms. Doing so requires 
access to the Hodge decomposition of /, but it is rather restrictive to assume access 
to this decomposition a priori. We are able to access the Hodge decomposition of / 
weakly in our estimators below, but at the expense of requiring more regularity of / 
than is needed to write the Hodge Laplace problem (cf. §4.2). In the a priori setting 
it is often possible to obtain improved error estimates by considering the discrete 
variables and measuring their error separately in weaker norms; cf. Theorem 3.11 
of [4] . This is an interesting direction for future research in the a posteriori setting 
as it may help to counteract this "Hodge imbalance" in the residual. 

The paper is organized as follows. In Section 2 we review the Hilbert complex 
structure employed in flnite element exterior calculus, begin to develop a posteri- 
ori error estimates using this structure, and establish a framework for bounding 
errors in approximating harmonic forms. In Section 3, we recall details about the 
de Rham complex and also prove some important auxiliary results concerning com- 
muting quasi-interpolants and regular decompositions. Section 4 contains the main 
theoretical results of the paper, which establish a posteriori upper bounds for er- 
rors in approximations to the Hodge Laplacian for the de Rham complex. Section 5 
contains corresponding elementwise efficiency results. In Section 6 we demonstrate 
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how our results apply to several specific examples from the three-dimensional de 
Rham complex and where appropriate compare our estimates to previous literature. 

2. HiLBERT COMPLEXES, HARMONIC FORMS, AND ABSTRACT ERROR ANALYSIS 

In this section we recall basic definitions and properties of Hilbert complexes, 
then begin to develop a framework for a posteriori error estimation. 

2.1. Hilbert complexes and the abstract Hodge Laplacian. The definitions 
in this section closely follow [4], which we refer the reader to for a more detailed 
presentation. We assume that there is a sequence of Hilbert spaces W'^ with inner 
products (•,•) and associated norms || • || and closed, densely defined linear maps 
d'' from W'^ into VF'^+i such that the range of d'^ lies in the domain of 6?*^+^ and 
^fe+i o d'^ = 0. These form a Hilbert complex {W,d). Letting V'' C W'^ be the 
domain of d'' , there is also an associated domain complex (V, d) having inner prod- 
uct {u,v)yk — {u,v)yyk + {d^u, d^ v)-^k+i and associated norm || • \\v- The complex 
... — > V'^~^ — 7- F*^ — > V'^^^ -^ ... is then bounded in the sense that d'^ is a bounded 
linear operator from V*' into V^'^^ . 

The kernel of d^ is denoted by 3*^ — S'' © ^^ ^ where S*^ is the range of d'^"^ 
and S)^ is the space of harmonic forms '<&^^^ n ^^ . The Hodge decomposition is an 
orthogonal decomposition of W^ into the range 55*^, harmonic forms S)^ , and their 
orthogonal complement 3*^^"' . Similarly, the Hodge decomposition of V^ is 

(2.1) V^ ^%''®S)''®-i^^, 

where henceforth we simply write 3*^^ instead of 3^^^^ except as noted. The dual 
complex consists of the same spaces W^ , but now with increasing indices, along 
with the differentials consisting of adjoints d*^. of d^~^ . The domain of d*^ is denoted 
by 14*, which is dense in W^ . 

The Poincare inequality also plays a fundamental role; it reads 

(2.2) ll«lk<l|rf't'lk, «e3'^. 

Here and in what follows, we write a < & when a < Cb with a constant C that 
does not depend on essential quantities. Finally, we assume throughout that the 
complex {W,d) satisfies the compactness property described in §3.1 of [4]. 

The immediate goal of the finite element exterior calculus framework presented in 
[4] is to solve the "abstract Hodge Laplacian" problem given by Lu — {dd* +d* d)u — 
f. L : W'' -^ W^ is called the Hodge Laplacian in the context of the de Rham 
complex (in geometry, this operator is often called Hodge-de Rham operator) . This 
problem is uniquely solvable up to harmonic forms when f J- Sj'^. It may be 
rewritten in a well-posed weak mixed formulation as follows. Given / G W'', we 
let p = Pfjkf be the harmonic portion of / and solve Lu ~ f ^ p. In order to 
ensure uniqueness, we require u _L Sj''. Writing a — d*u, we thus seek {a,u,p) € 
V''^^ xV'' X S}'' solving 

(cr, r) - (dr, m) === 0, TeV''-\ 

(2.3) {dc7,v) + {du,dv) + {v,p) == {f,v), l;ey^ 

{u,q) = 0, qe^\ 

So-called inf-sup conditions play an essential role in analysis of mixed formula- 
tions. We define i?(i7,u,p;r,w, g) — {a,T) — {dT,u) + {da,v) + {du,dv) + {v,p) — {u,q), 
which is a bounded bilinear form on [V''^^ xV'' x Sj''] x [V''~^ x V'^ x Sj'^]. We will 
employ the following, which is Theorem 3.1 of [4]. 
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Theorem 1. Let {W,d) be a closed Hilbert complex with domain complex [V^d). 
There exists a constant 7 > 0, depending only on the constant in the Poincare 
inequality (2.2), such that for any [a^u^p) G V^^^ x V^ x S^^ there exists (T,v,q) G 
l/*^-! xV'^ X Sj'^ such that 

(2.4) B{a,u,p;T,v,q) > 7(|k||y + \\u\\v + b||)(!|T||v + \\v\\v + \\q\\). 

2.2. Approximation of solutions to the abstract Hodge Laplacian. Assum- 
ing that {W, d) is a Hilbert complex with domain complex (V, d) as above, we now 
choose a finite dimensional subspace V^f C V'^ for each k. We assume also that 
dV^ C V^^'^ , so that {V/^,d) is a Hilbert complex in its own right and a subcom- 
plex of (y, d). It is important to note that while the restriction of d to V^ acts as 
the differential for the subcomplex, d* and the adjoint d*f^ of d restricted to V^ do 
not coincide. The discrete adjoint djj does not itself play a substantial role in our 
analysis, but the fact that it does not coincide with d* should be kept in mind. 
The Hodge decomposition of V^ is written 

(2.5) V/:^^t®Sj1^®3'H^. 

Here *8^ = dV/^, with similar definitions of Jo^ and 3^^ where _L is in Vh- This 
discrete Hodge decomposition plays a fundamental role in numerical methods, but 
it only partially respects the continuous Hodge decomposition (2.1). In particular, 
we have: 

(2.6) iof; C 3' but io^ ^ ^^ 

3^ ^ 5"^- 



Bounded cochain projections play an essential role in finite element exterior 

h 



calculus. We assume the existence of an operator tt^ : V*^ — )■ Vt which is bounded 



in both the ly-norm || • || and the V^-norm || • ||v and which commutes with the 
differential: d^ o tt,^ ~ ''^h^^ o d'^ . In contrast to the a priori analysis of [4], our a 
posteriori analysis does not require that Tih be a projection, that is, we do not require 
that -Kh act as the identity on Vt- In more concrete situations we shall however 
require certain other properties that are not needed in a priori error analysis. 

Approximations to solutions to (2.3) are constructed as follows. Let {uh^ '^hiPh) G 
Vl^-^ X Vh X ^h satisfy 



(2.7) {dah,Vh) 



Existence and uniqueness of solutions to this problem are guaranteed by our as- 
sumptions. A discrete inf-sup condition analogous to (2.4) with constant 7^, depend- 
ing on stability constants of the projection operator tt/j but otherwise independent 
of Vh is contained in [4]; we do not state it as we do not need it for our analy- 
sis. In addition. Theorem 3.9 of [4] contains abstract error bounds: So long as the 



{crh,Th) - {dTh,Uh) = 


= 0, 




{duh,dvh) + {vh,Ph) = 


= (/,«/.), 


Vh e VI:, 


{uh.qh) = 


= 0, 


Qh e Sjl 
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subcomplex {Vh,d) admits uniformly V^-boundcd cochain projections, 

\\o- - o-hWv + \\u- Uh\\v + \\p-Ph\\ 

< inf llcr — t||i/+ inf IIm — flli/ 
(2.8) rev^' -ev/; 

+ inf \\p - q\\v + fj- inf \\P<bu- v\\v, 

where jl — sup^^^^fc ||r||=i II (-^^""ft)^!!- We will use the notation Ps for the orthogonal 
projection onto the subspace S as in the case of Peg above. 

2.3. Abstract a posteriori error analysis. We next begin an a posteriori error 
analysis, remaining for the time being within the framework of Hilbert complexes. 
A working principle of a posteriori error analysis is that if a corresponding a priori 
error analysis employs a given tool, one looks for an a posteriori "dual" of that 
tool in order to prove corresponding a posteriori results. Thus while the proof of 

(2.8) employs a discrete inf-sup condition, we shall employ the continuous inf-sup 
condition (2.4). Writing e^ = a — a^,, e„ = m — u^, and Cp = p — ph, we use the 
triangle inequality and (2.4) to compute 

||eo-||y + ||e„||v + ||ep|| < (||eo-||y + ||e„||y + \\p - PsjPhW) + \\P?,Ph ^ Ph\\ 

<- sup B(e^,eu,p- P<r,Ph;T,v,q) + \\P^ph-ph\\ 

^ (r,t),9)ey'=-ixy'=xi5'', 

|T||v + ||f II V + lklHl 

(2.9) 1 / 

< - sup ( (e^, t) - (dr, e„) + {de„, v) + (de„, dv) 

'y {T,v,q)eV''-'^xV''xS^'', ^ 

\\T\\v + \\v\\v + \\q\\ = i 



+ {v,ep) + {eu,q)) + (1+ -)\\P^Ph -Ph\ 



7 

Employing Galerkin orthogonality implied by subtracting the first two lines of (2.7) 
and (2.3) in order to insert tt^t and tt/jW into (2.9) and then again employing (2.3) 
finally yields 

||ecr||v+||e„||y + j|ep|| 

<- sup Uah.T -TThT) - {d{T-TThT),Uh) 

7 (t,v ,q)eV^-^ y.V'' y.^'' , ^ 

|t|| v + ||'u||v + ||<j||=l 



(2.10) 



+ (./ - dcFh -ph,v- TThv) - {duh, d{v - TThv)) + (e„, q)j 



+ il + -)\\Ps,Ph-Ph\\- 

7 

The terms {ah,T-TThT) -{d{T-TrhT),Uh) and {f -d(Th-Ph,v-'Khv) + {duh,d{v- 
TT/jw)) in (2.9) can be attacked in concrete situations with adaptations of standard 
techniques for residual-type a posteriori error analysis, but no further progress can 
be made on this abstract level without further assumptions on the finite element 
spaces. The terms (e„,(7) and (1 + -)\\Pf^ph — Ph\\j on the other hand, are nonzero 
only when i^^ ^ Sj'^. In this case (2.7) is a generalized Galerkin method, and further 
abstract analysis is helpful in elucidating how these nonconformity errors may be 
bounded. We carry out this analysis in the following subsection. 
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2.4. Bounding the "harmonic errors". We next lay groundwork for bounding 
the terms \\ph-PsoPh\\ and supgg^fc(e„, g). Since p^ ^ ^h^ (2-6) and S'' ==<B*=®.^'' 
imply that p^ — P^ph € ^^ ■ Recalling that v e ^^ implies that v — dcf) for some 
(j) e V^~^ and also that PfyPh e S^^ ^-^^ yields 

(2.11) \\pf^-P^p^\\^ sup (p^- Pf^p^^v) ^ sup {ph,d(j>). 

t>e®'=,||«||=i 0ey=-i,||d0||=i 

The discrete Hodge decomposition (2.5) implies that n^dcl) — dir^^^^cf) G 25^ _L i^Jj 9 
Ph- Also note that by the Poincare inequality (2.2), s.\xp^^v''-^.\\d(t>\\=i{Phid4') is 
uniformly equivalent to sup^gyt-i ii^ii ^^^(p/i, di/)). Thus 

(2.12) \\Ph-P^Ph\\< sup {ph,d{(j)-'Kh(j})). 

0ev'=-MI0||v=i 

We do not manipulate (2.12) any further without making more precise assumptions 
about the spaces and exterior derivative involved. Recall that the goal of (2.12) 
is to measure the amount by which the discrete harmonic function p^ fails to be a 
continuous harmonic function. If p^ were in fact in Jo'', we would have d*ph = 0, 
which would immediately imply that the right-hand-side of (2.12) is 0. In (2.12) we 
measure the degree by which this is not true by testing weakly with a test function 
d^, minus a discrete approximation to the test function. 

Before bounding the term sup g^h \\q\\=i{GmQ) we consider the gap between S)^ 
and S)^. Given closed subspaces A, _B of a Hilbert space W , let 

(2.13) 5{A,B)^ sup dist(a;,B)= sup ||x-PBa;|l. 

a:G A,||a:|| — 1 2;GA, ||a;|| — 1 

The gap between the subspaces A and B is defined as 

(2.14) gap(A, B) = max((5(A, B),S{B, A)). 

In the situation below, we will require information about (5(io'^,ioJj), but are able 
to directly derive a posteriori bounds only for S{S)'j^,S)''). Thus it is necessary to 
understand the relationship between S{A,B) and S{B,A). 

Lemma 2. Assume that A and B are subspaces of the Hilbert space W , both having 
dimension n < oo. Then 

(2.15) 5{A,B)^5{B,A) = g&^{A,B). 

Proof. The result is essentially found in [18], Theorem 6.34, pp. 56-57. Assume 
first that 5{A,B) < 1. The assumption that dim A = dimS then implies that 
the nuUspace of Pb is and that Pb maps A onto B bijectively. Thus Case i of 
Theorem 6.34 of [18] holds, and (2.15) follows from (6.51) of that theorem by noting 
that 5{A,B) = \\I-Pb\\{a,w) = ||(/ - Pb)^^!!^^,^). 

If 5{A,B) ~ 1, then there is 7^ fe G S which is orthogonal to Pb{A). Let- 
ting {ai}i=i,...M be an orthonormal basis for A, wc have P46 = X]i=i('^i' ^)'^i ~ 
Tf=l{PBa^,b)a, = 0, since 6± Pb(A). Thus 1 = \\I - Pa\\(b,A) = KB, A), so that 
(2.15) holds in this case also. D 

Thus we can bound gdsp{S)^ ^ S)\) by bounding only 5{S)\,S)^), which we now turn 
our attention to. First write 5{S)^,S)^) = sup^^g^^fe 11^^11^]^ \Wh — P^kqh\\. For a given 
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qh e -QJJ, we may employ exactly the same arguments as in (2.11) and (2.12) above 
to find 

(2.16) \\qh-P9}qh\\< sup {qh,d{(p-TTh(f>)). 

We now let {qi, ..., qm} be an orthonormal basis for Jojj and assume that we have a 
posteriori bounds 

(2.17) sup (g„d((/)-7r,,0))<^„ z = l,...,M. 

0ey'=-i,||0||v=i 

We obtain such bounds for the dc Rham complex below. Given an arbitrary unit 
vector qfi G S)^, we write qh — X]i=i '^iQi^ where \a\ = 1. Inserting this relationship 
into (2.16) yields S{Sj'j^,Sj'^) < supig^i j^ X]i=i ^^il^i- This expression is maximized by 
choosing a = /i/|/i|, where fl ~ {/ii, .., /im}- Thus 

(2.18) S{^l^'^)<\^I\. 
Combining (2.18) with (2.15), we thus also have 

(2.19) gap(i3^i^;;)<|A^|. 

Now we turn our attention to bounding ||Pf|U;i|| ~ sup g^t j|„j|=i(e„, g). Our 
analysis of this term is slightly unusual in that we suggest two possible approaches. 
One is likely to be sufficient for most applications and is less computationally in- 
tensive. The other more accurately reflects the actual size of the term at hand, but 
requires additional computational expense with possibly little practical payoff. 

We first describe the cruder approach. Because u^ _L i^^, 

\\Ps^Uh\\ ^ sup {q^uh) ^ sup {q- Ps^kq,Uh) 
(2 20) <?e-*i'°JI<?ll=i 'je-f3'=,||g||=i 

<<5(i3^i^|;)h,lHgap(^^i3^)h,lj. 

(2.19) may then be used in order to bound ga,p {Sj'^ , Sj'j^) . 

Next we describe the sharper approach. Since u^ _L Sj^, we have u^ = u^ + uj^, 
where u^ € *B^ and u^ £ 3^^- Since *B^ C S'' -L fj'', m^ -L Sjf^, and i^^ and S)'' 
are both perpendicular to 3*^^, we thus have for any q E Sj'' with ||g|| = 1 that 

, {uh,q) = {ui,q) = {ui,q-Ps^,ji) = {u^ - P^±uj^ ,q - P^^q) 

< \\ui - P^.^uiWWq - Ps,,q\\ < gM^',&^i - ^3-411- 

But 

(2.22) M^ - P^^ui = P^ui + P^ui = P^ui + P^uh. 

Here the relationship Ps^uj^ = PsjUh holds because *B^ C 25*^ and so Ps^Uh — 0. 
Thus ||m^ - Ps^M^li < \\P<sui:\\ + WPs^UhW- ||-P<8W^|i may be bounded as in (2.12) 
and (2.16) above: 

(2.23) 



\P<sUh\\= sup (u,, , d(</) - 7r/i(/))) 

< sup {uj^,d((l)-TTh(p)). 

0ev"=-i,||</.||v=i 



Assuming a posteriori bounds ga.p{^'^ , Sjf^) < fi and ||P<sm^|| < e, we thus have 

(2.24) \\P,,Uh\\<fi{^+\\Pf,Uh\\). 
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Inserting (2.20) into (2.24) then finally yields 

(2.25) \\P,,Uh\\<e|l + ^l''\\uh\\. 

We now discuss the relative advantages of (2.20) and (2.25). The corresponding 
term in the a priori bound (2.8) is /iinf^gyfc ||^®u ~ ^llvi which is a bound for 
||P(j^u|| (note the symmetry between the a priori and a posteriori bounds). The 
term jl (defined following (2.8)) is easily seen to be bounded by gap(i3'^,io{5) at 
least in the case that -Kh is a VF-bounded cochain projection. Also, it is easily seen 
that jl is generally of the same or higher order than other terms in (2.8) when 
standard polynomial approximation spaces are used. Carrying this over to the a 
posteriori context, (2.20) will yield a bound for HPj^u^H that while crude is not 
likely to dominate the estimator or drive adaptivity in generic situations. 

If a sharper bound for UPsjU/iH proves desirable (e.g., if gap(Jo'',.^^J||u/i|| dom- 
inates the overall error estimator), then one can instead employ (2.25). This cor- 
responds in the a priori setting to employing the term inf^gyk Ij^fs^ ^ ^||y and 
is likely to lead to an asymptotically much smaller estimate for |JPqM/i||. How- 
ever, computing the term e in (2.25) requires computation of the discrete Hodge 
decomposition of Uh, which may add significant computational expense. 

2.5. Summary of abstract bounds. We summarize our results above in the 
following lemma. 

Lemma 3. Assume that (W,d) is a Hilbert complex with suhcom,plex (V/i,d) and 
commuting, V -bounded cochain operator nh '■ V — ?• Vh, and in addition that {a,u,p) 
and {ah,Uh,Ph) solve (2.3) and (2.7), respectively. Then for some {T,v,q) £ V'^^^x 
V^ X 9)^ with \\t\\v + \\v\\v + ll^ll = 1 and some cj) <E V^^^ with \\(i>\\v ~ 1, 

116^11^ + ||e„||v + ||ep|| < |(e^,T-7r/jr) - (d(T - tt/it), e„)| 

(2.26) + l(./ - d<^h -Ph,v- TThv) - {duh,d{v - tt^w))! 
+ \{Ph,d((j}-TTh4>))\+ ^l\\ui -PjiU^II. 

Here ^ ^ (Ej=iMf)^^^; where sup^^yk-i^\\^^^i{qi,d((f> ~ iThcf))) < Mj for an or- 
thonormal basis {qi, ....,gAf} of S)^. For the last term in (2.26) we may either use 
the simple bound \\uj^ — P^±uj^\\ < \\uh\\ or employ the bound /x||u^ — P-^±uj;^\\ < 
/ie + /i^ 1 1 u^ 1 1 , where 

(2.27) sup {ui,d{^-TTh^))<e. 

3. The de Rham complex and commuting quasi-interpolants 

As above, we for the most part follow [4] in our notation. Also as above, we shall 
often be brief in our description of concepts contained in [4] and refer the reader to 
§4 and §6 of that work for more detail. 

3.1. The de Rham complex. Let J7 be a bounded Lipschitz polyhedral domain 
in M", n > 2. Let A'^(i7) represent the space of smooth /c-forms on il. A'^(f7) is 
endowed with a natural L2 inner product (•, •) and L2 norm || • || with corresponding 
space L2A'^(J7). Letting also d be the exterior derivative, HA'^{fl) is then the 
domain of d'' consisting of L2 forms fi for which doj G _L2 A'^+^ (fi) ; we denote by || • || /f 
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the associated graph norm. {L2A {n),d) forms a Hilbert complex (corresponding 
to (W, d) in the abstract framework of the preceding section) with domain complex 

(3.1) 0^i7A°(O) A HK^{9.) -A---A i7A"(f7) ^ 

corresponding to (V, d) above. In addition, we denote by VFIA'^(i7) the correspond- 
ing Sobolev spaces of forms and set H^A''{n) = VFJA'^(ri). Finally, for cj C M", we 

let II • \\lo = II • ||L2A'=(aj) and || • \\h,ui = \\ ■ \\hai'{lj)', in both cases we omit a; when 
w = il. 

Given a mapping (p : fti ^ ^2, we denote by (p*uj G A'^(r2i) the pullback of 
uj e A''(rJ2), i.e., 

(3.2) {(l)*Oj)x{vi,...,Vk) = UJ^(^^){D(l)^{vi), ..., D(pa:{vk)). 

The trace tr is the pullback of w from A''{fl) to A'^{dn) under the inclusion, tr is 
bounded as an operator HA''{n) ~> H-^/^A''{dn) and H^A''{n) -^ H^^^A''{dn), 
and thus also H^A''{n) -^ L2A^{dVL). We may now define HA''{n) = {cj e 
HA^{ijj) : trw = on 9il}. In addition, we define the space HlA*'{Vt) as the 
closure of C^A^{Vt) in H^A^{Vt). Hq A'^(f2) essentially consists of forms which are 
in every component on 90, which is in general a stricter condition than trw = 0. 
The wedge product is denoted by A. The Hodge star operator is denoted by * 
and for u; e A*^, ^ e ^y""*^ satisfies 



(3.3) wAm= (*^,m)vo1, / wAm= (*'^,m)l2A"-'=(!:^o)- 

-k is thus an isometry between L2A^ and L2A"^'^. The coderivative operator S : 
j^k _^ pk-i jg defined by 

(3.4) i<5uo = {-lfd-kuj. 

Applying Stokes' theorem leads to the integration-by-parts formula 

(3.5) (duj,^) ^ (uj.Sfi) + trwA tr Tk-^, w e i/A''"^ ^ e iJ^A''. 

Jan 

The coderivative coincides with the abstract codifferential introduced in §2.1 when 
tr QQ -k ^ = 0. That is, the domain of the adjoint d* of d is the space 7J*A'^(il) 
consisting of forms /i e L2A'' whose weak coderivative is in L2A''^^ and for which 
tr • /J. = 0. We will also use the space H* A^ — -k{HK"-^^) consisting of L2 forms 
whose weak codifferential lies in L2; note that v G H* A'^ implies that tr -kv e H^^/^ . 
The Hodge decomposition L2A^{Vt) == ^B*^ Jo'' ®^ consists of the range S'^ = 
{dip : ip e HA^~^{Vl)}, harmonic forms S^^ = {uj ^ HA^{Vt) : dw = 0, 5uj = 
0, tr *a; = 0}, and range <B^ = {(5w : w G i?*A'=+i(fi)} of 5. dimi^*^ is the k-ih. 
Betti number of fJ. The mixed Hodge Laplacian problem corresponding to (2.3) 
now reads: Find (cr, M,p) G HA'^^^ x HA'' x S)'' satisfying 

(3.6) a — du, da + Sdu = / — p in f2, 

(3.7) tr*u = 0, tr *du^ on dn, 

(3.8) u ± ^^ 

The boundary conditions (3.7) are enforced naturally in the weak formulation (2.3) 
and so do not need to be built into the function spaces for the variational form. 
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The additional boundary conditions 

(3.9) tr • cr = 0, ti ir5du = on dfl 

are also satisfied. To see this, note that d and tr commute since tr is a puUback, 
and that tr -ka and tr -kSdu are both well defined in H^^^^ since 6a = and SSdu = 
imply that (T,Sdu e H* . Thus by (3.4), tr -ka = tr {-l)''d*u = (-l)''d tr -ku = 0. 
Similarly, tr * 6du = tr (— 1)*^^*^^ — (— l)'^(itr -kdu = 0. These relationships are 
roughly akin to noting that for a scalar function m, the boundary condition u = 
on dn implies that the tangential derivatives of u along dfl are also 0. 

We also consider the Hodge Laplacian with essential boundary conditions, that 
is: Find {a,u,p) G HK^^^ x Hk^ x S)^ satisfying 

(3.10) a — 5u, da + Sdu ^ f — p in fl, 

(3.11) tro- = 0, tr u = on 9f7, 

(3.12) u _L ^^ 

This is the Hodge Laplace problem for the de Rham sequence (3.1) with each 
instance of HA''{Q) replaced by HA'^{V,). Here we have denoted the correspond- 
ing parts of the Hodge decomposition using similar notation, e.g., Sj'' = {oj € 
3'=|(w,Ai) -0, M £»'=}. 

3.2. Finite element approximation of the de Rham complex. Let Th be a 

shape-regular simplicial decomposition of O. That is, for any Ki,K2 € Th, Kir]K2 
is either empty or a complete subsimplex (edge, face, vertex, etc.) of both Ki and 
K2, and in addition all K £ Th contain and are contained in spheres uniformly 
equivalent to hx '■— diani(ii'). 

Denote by (V/i, d) any of the complexes of finite element differential forms con- 
sisting of Vr and V~ spaces described in §5 of [4]. We do not give a more precise 
definition as we only use properties of these spaces which are shared by all of them. 
The finite element approximation to the mixed solution (cr, u, p) of the Hodge Lapla- 
cian problem is denoted by (ah, Uh,Ph) & V^ x V^~^ x^^ and is taken to solve (2.7), 
but now within the context of finite element approximation of the de Rham com- 
plex. In order to solve (3.10)-(3.12) we naturally employ spaces V/^ — V/^ O H A'^ (i^) 
of finite element differential forms. 

3.3. Regular decompositions and commuting quasi-interpolants. We also 
employ a regular decomposition of the form lu = dip + z, where lo G HA only, but 
93, z e H^. In the context of Maxwell's equations the term "regular decomposition" 
first appeared in the numerical analysis literature in the survey [17] by Hiptmair 
(although similar results were previously available in the analysis literature) . Pub- 
lished at about the same time, the paper [22] of Pasciak and Zhao contains a similar 
result for if (curl) spaces which is also often cited in this context. Below we rely 
on the paper [20] of Mitrea, Mitrea, and Monniaux, which contains regularity re- 
sults for certain boundary value problems for differential forms that may easily be 
translated into regular decomposition statements. Recent work of Costabel and 
Mcintosh [14] contains similar results for forms, though with handling of boundary 
conditions that seems slightly less convenient for our purposes. 

We first state a lemma concerning the bounded invertibility of d; this is a special 
case of Theorem 1.5 of [20]. 
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Lemma 4. Assume that B is a bounded Lipschitz domain in M" that is homeomor- 
phic to a ball. Then the boundary value problem dtp = g E -L2A'^(_B) in B, tr (^ = 
on dB has a solution (f G HqK^~^{B) with \\^p\\h'^k''-'^{b) ^ llffll-B */ o,nd only if 
dg = in B, and in addition, tr g — on dB ifO<k<n— 1 and jgg = ^ if 
k = n. 

Employing Lemma 4, we obtain the following regular decomposition result. 

Lemma 5. Assume that CI is a bounded Lipschitz domain in M", and let < k < 
n - 1. Given v e Hk^{VL), there exist Lp G H^K^^'^(Vl) and z G H^A''{^) such that 
V — dip + z, and 

(3-13) ||</'||HiAfc-i(o) + \\z\\H^A''{n) ^ I|w||ha'=(o)- 

Similarly, if v e HA''{i}), then there exist ip <E H^A'^^^iil) and z G _ffgA'^(ri) such 
that V = dip + z and (3.13) holds. In the case k ^ n, HA"'{Q,) is identified with 
L2A"(r2). The same results as above hold with the exception that z G Lit^{^ only 
and satisfies H^IU^A^cn) ^ hWHA^iQ)- 

Proof. We first consider the case v G HA''{il). By Theorem A of [21], the assump- 
tion that dfl is Lipschitz implies the existence of a bounded extension operator 
E : HA'' {ft) -^ HA''{R"). Without loss of generality, we may take Euj to have 
compact support in a ball B compactly containing fl, since if not we may multiply 
Ev by a fixed smooth cutoff function that is 1 on 51 and still thus obtain an HA- 
bounded extension operator. Assuming that < k < n — 1, we solve dz = dEv for 
z G Hq{B) and dip = Ev — z ior ip £ HqA''^^{B), as in Lemma 4. The necessary 
compatibility conditions may be easily verified using rf o d = 0, rftr = trd, and 
dEv = dz. Restricting ip and z to 51, we obtain (3.13) by employing the bounded- 
ness of E along with Lemma 4. In the case k = n, we let z = {\B\~^ J^ Ev)vol, 
where vol is the volume form. We then have Jg{Ev — z) = 0, and proceeding by 
solving dip — Ev — z as above completes the proof in this case also. 

In the case v G i?A'^(il), Lemma 5 may be obtained directly from Lemma 4 when 
n is simply connected by applying the procedure in the previous paragraph with 
B = Q. The general case follows by a covering argument. Let {^i} be a finite open 
covering of f2 such that 11 n 51^ is Lipschitz for each i, and let {xi} be a partition 
of unity subordinate to {51^}. When < A: < n — 1, we first solve dzi — d{xiv) for 
Zi G iJoA'^(ili n ft) and let z ~ J^n ^i ^ ^^^''{Vt). A simple calculation shows 
that the compatibility conditions of Lemma 4 are satisfied, and in addition dz ~ dv 
since X^n Xi ~ 1- Similarly, we solve dpi ~ Xii"^ ^ ^) ^oi ipi G HQA''^^{ili) and 

set ip — X^o fi- -'-^ ^^^ ''9'S^ k ^ n, let Zi = I |Hi n il|^^ In nnXi'^) "^ol and let 
ipi G HQA"^^{fl n fli) solve dipi = XiV — Zi. Setting p — X^n Vi^ z — v — dip, and 
employing Lemma 4 completes the proof. D 

Our next lemma combines the regular decomposition result of Lemma 5 with a 
commuting quasi-interpolant in order to obtain approximation results suitable for 
a posteriori error estimation. Schoberl defined such an interpolant in [23] for the 
classical three-dimensional de Rham complex and extended his results to include 
essential boundary conditions in [24]. We generally follow Schoberl's construction 
here, although our notation appears quite different since we use the unified notation 
of differential forms. Schoberl instead employed classical notation, which makes 
the necessary patterns clear and perhaps more concrete but also requires a different 
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definition of the intcrpolant for each degree of forms (H^, H{cur\), iJ(div), and L2). 
Our development also has many similarities to that of Christiansen and Winther 
in [13], who develop a commuting projection operator for differential forms. Their 
operator is however global and not suitable for use in a posteriori error estimation. 
In [24] Schoberl defines and analyzes a regular decomposition of the difference 
between a test function and its interpolant over local element patches for the three- 
dimensional de Rham complex, whereas we first carry out a global regular de- 
composition and then interpolate. The only advantage of Schoberl's approach in 
the context of a posteriori error estimation seems to be localization of domain- 
dependent constants, which appears to be mainly a conceptual advantage since the 
constants are not known. Also, commutativity of the quasi-intcrpolant is not nec- 
essary for the proof of a posteriori error estimates, although it may simplify certain 
arguments. See [12] for proofs of a posteriori estimates for Maxwell's equations that 
employ a global regular decomposition but not a commuting interpolant. Commu- 
tativity may however be useful in other contexts. For example, one may modify 
the proofs of quasi-optimality of an AFEM for Maxwell's equations in [27] to use 
a global regular decomposition and a quasi-intcrpolant instead of a local regular 
decomposition of the interpolation error so long as the interpolant commutes. 

Lemma 6. Assume thatv G _ffA'^(il) with \\v\\h < 1- Then there exists an operator 
nf' : i2A*'(ri) -> V^ such that d'^+^II^ = n|;+^d*', and in addition the following 
hold. If k — 0, HA'^ = H^ holds and we have 

(3.14) J2 ['^K'h-^hv\\l + h],''\\tTiv-Ilhv)\\li, + \v-nhv\l,f^j,-^'\<l. 
Ken 

If^<k< n-1, there exist (p G H^A''^'^{n) and z e i/^A*'(ri) such thatv = dcp+z, 
Uf^v = dnj;" V + ^h^: and for any K £%, 

J2 [h],'{y-nhif\\l + \\z-nhz\\l) 

(3.15) ^67; 



/i^^dl tr (^ - n;,^)||2^ + II tr (z - n;,z)||2^) 



< 1. 



In the case k — n, the space HA''{il) is L2A^{Vl), and there exist Lp E H^A^ ^{^)j 
z G i2A"(il) such that v ~ dip + z, Il^f = dWhip + II^z, and 

HK) + ll^ ^ n/iz||^^^fc(^)) 
(3.16) 

' < 1- 

Assume that 1 < k < n and cj) G HA''~^{fl) with \\(J>\\h < 1- Then there exists 
p G -ff^A'^~^(ri) such that dp — d(j>, Hhdcj) = dllh(f> = dllfi(p, and 



<en 


k'QI^ 


-'^iMll^^A''-^ 


+h,'\ 


tT{ip- 


-^hp)\\l^(OK) 



Y. [''kW^ - ^^'^Wk + v'li tr {p - ^hpWdK 

Ken 



< 1. 



Finally, the above statements hold with IIA''{n) replaced by IlQA'^{n) and _ff^A*'(ri) 
replaced by HqA'^^U), and in this case Uh : HA'^ifl) -^ Vj^ . 

Proof. Let lih = IhRhi where following [24, 13] Ih is the canonical interpolant for 
smooth forms and R\ is a smoothing operator with smoothing parameter e which 
we detail below. (Note that the construction in [13] involves a further operator J^, 
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which is the inverse of Hh apphcd to V/^ . J^ is non-local and thus not suitable for 
use here.) We conceptually follow [24] in our definition of i?| but also use some 
technical tools from and more closely follow the notation of [13]. We omit a number 
of details, so a familiarity with these works would be helpful to the reader. 

Following [13], there is a Lipschitz-continuous vector field X{x) defined on a 
neighborhood of Q such that X{x) ■ n{x) > for all outward unit normals n{x), x € 
dfl. Let gh{x) be the natural Lipschitz-continuous mesh size function. There is then 
(5 > so that for e > sufficiently small depending on 51, B^{x + 6egh{x)X(x)) C 
M"\ri for all x e n with dist(x,9r2) < egh{x), and Bf{x - Segh{x)X {xj) C ft for all 
X G dil. Let xv be the standard piecewise linear "hat" function associated to the 
vertex V, and let A/g be the set of boundary vertices. Extending xv and g^ by re- 
fiection over 357 via X (cf. [13]), we define ^d(x) ~ ^ + ^veJ\f Xv {x)Segfi(x) X (x) 
and <^n{x) =x- EyeTVg Xv{x)6egh{x)X{x). 

Following [24], we associate to each vertex V in Th a. control volume ily- Let fly 
be the ball of radius egh{V)/2 centered at 1/ if t^ is an interior or Dirichlet boundary 
node, and at ^n{V) if y is a boundary Neumann node. Also let fv G P„ satisfy 



(3.18) 



p{x)fv{x)dx =p{V), pe 




Figure 1 . Schematic of transformations and control volumes for 
the Dirichlet (left) and Neumann (right) cases. 



Given K G Th, ict {Vi, ...,Vn+i} be the vertices with associated barycentric co- 
ordinates {Ai(x),...,A„+i(a;)} for x € K. Let f (x, j/i, ....,i/„+i) = Ya^^ Xi{x)yi = 
X + YL^=i ^i{x){,yi — Vi) for Hi € Oy. . We also denote by w the extension of w by 
to M" in the case of Dirichlet boundary conditions, and the extension of w to a 
neighborhood of J7 by smooth reflection via X in the case of Neumann boundary 
conditions (cf. [13]). Letting 2: = x in the case of natural boundary conditions and 
X — $£) o i; in the case of essential boundary conditions, we define 



(3.19) {Rlu:), 



Vl 



fvi ivi) ■ ■ ■ /y„+i {V7i+i){x*uj)^ dy„+i • • • dyi. 



Commutativity of i?^ and thus of li^ immediately follows as in [23, 24, 13]. In 
addition, w = on M" \ Vt implies that {Rj^uj)^ = for a; G dil in the case of 
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essential boundary conditions, since in this case d[st{x{x),dfl) < gh{x)e for all 

(2/l,---,l/n+l) e rJvj X ••• X ^v„^^_ SO that $£,of(a;,yi,...,y„+i) G M"\f7. 

We now establish that i?^ preserves constants locally on elements sufficiently 
removed from 90. This implies the same of XI/i since /^ preserves constants. Let 
w be a constant fc-form, and let a; G fJ lie in K Cz Th which in the Neumann 
case satisfies K n Sfi = and which in the Dirichlet case satisfies ujk H dfi = 0. 
Noting that now x = x on ujk, computing that Dx{x, j/i, ..., yn+i) = I + Y!^^=i iVi ^ 
Vi) <Si VAi(x), applying (3.2), and applying the multilinearity of cj, wc find that for 
n- vectors wi, ..., w^ 

{Rh^)x{vi,...,Vk)= I ••/ .M(yi)---/y„+i(2/n+i) 

/ n+l n+1 \ 

(3.20) X w Li + ^(VA,(a;) • vi){y, - K), ...,Wfe + ^(VA,(x) • Vk){y^ - Vi) 

\ i=l i=l J 

X dyn+i • • • dj/i 

= Uj{vi,...,Vk) +$. 

Here $ consists of a sum of constants multiplying terms of the form lu(zi, ..., Zk), 
where Zi = yj — Vj for at least one entry vector Zi and for some 1 < j < n + 
1. (3.18) then implies that /j-^^^ • • • /j-^^^ ^ fv^ (yi) ■ ■ ■ fv„+i (j/n+i)* dy„+i • • • dyi = 

0. Observing that J^ fv{y)'^y ^ 1 completes the proof that R\uj{vi, ...,Vk) — 
Lo{vi, ..., Wfe) in the case that x lies in an interior element. 

Arguments as in [13, 23, 24] show that for a; G i2A'^, ||nft,c<;||/i- < ||a;||x*, where 
K* = ijJk if K is an interior clement and K* — lox U {x e M" : dist(a;, iiT) < ehx} 
otherwise; cf. Figure 4.2 of [13]. In the latter case the definition of the extension 
of a; to M" \ 51 as either or the puUback of uj under a smooth refiection implies 
that the values of lj on K* in fact depend only on values of w on ujk so long as e is 
sufficiently small, which in turn implies that ||n;iaj||/i- < jli^ll/f ^ ll'^llw^- For ii' not 
abutting 957 we may combine the above properties of n^j with the Bramble-Hilbert 
Lemma (cf. [9]) in order to yield h'^\\z -Iihz\\K + \z -Ii-hz\H^\''(K) ^ klHiA'=(w;c) 
for z G H^A^{Vl). A standard scaled trace inequality for z E H^K^{K) reads 

Ijtr (z - Tlhz)\\dK < h'^ \\z - n^zll^ + /i]/^|z - n/iz|^iAfc(K)- Combining these 
inequalities with the finite overlap of the patches ujk implies (3.15), (3.14), and 
(3.16) modulo boundary elements. 

When V^ = HA''{fl) and x e K with K n dil ^ ^, the above argument that 
constants are locally preserved and a standard Bramble-Hilbert argument applies 
holds so long as the convex hull of ilvi , ■■■,^v^+i lies in fl for the vertices Vi, .., Ki+i 
of K. This is true if hx < ho for some ho depending on fl. To prove this, let 
Vi, ...,Vk be the vertices of K lying on dfl. Our assumptions imply that Oy. C O, 
i = 1, ..., k. We must show that the convex hull of ^Vi, ■■■, ^Vk ^tlso is a subset of 
il. This set is the union of all balls -B£g^(j,)/2(y), where y = '^i^i ^i^ n (Vi) for 

< Ai < 1 satisfying J2i=i ^i — 1- Fixing such a y, let y — J2i=i ^i^i ^ dK n dfl. 
Then dist($jv(y),ari) > egh{y). But |$Ar(y) - y| < e5,,(2/)/iK||£'AC||L^, so that 
Begn(v)/2{y) C 57 SO long as /ik||£'^||l^(k) < 5. We thus take ho - ^ 



2||DX|| 



-Am 



Thus the results of Lemma 6 follow as above when hx < ho- If the convex hull of 
flviy •••7^y„+i does not lie in 57, then the integral (3.19) may sample values of 



ui^ 



for X ^ 57. Extension by pullback of a reflection does not preserve constant forms. 
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and so in this case II/i docs not necessarily preserve constants locally. The operator 
n^ is however still locally L2 bounded, and since in this case we have h^ > ft-o we 
may think of h^ as merely a constant depending on fl and thus still obtain the 
approximation estimates of Lemma 6 after elementary manipulations. 

In the case where V'' = HA''{Q) and x lies in a patch ojk with K n dV, ^ 0, 
Till does not preserve constants locally, but is still locally L2 bounded. Because we 
in this case apply Hh to forms z,ip E Hq (that is, to forms which are identically 
on dfl), we may apply a Poincarc inequality in conjunction with elementary 
manipulations in order to obtain (3.15), (3.14), and (3.16) in this case as well. 

(3.17) follows by a similar argument, that is, by extending (j) _ff-continuously to 
a ball B in the Neumann case, solving the boundary value problem dip — dcf) on 
B or fl as appropriate, employing the H^ regularity result of Lemma 4, and then 
applying properties of Hh- D 



4. Reliability of a posteriori error estimators 

In this section we define and prove the reliability of a posteriori estimators. We 
will establish a series of lemmas bounding in turn each of the terms in (2.26). Below 
we denote by |x] the jump in a quantity x across an element face e. In case e C dQ, 
Ixl is simply interpreted as x- AH of our results and discussion below are stated for 
the case of natural boundary conditions; results for essential boundary conditions 
are the same with the modification that edge jump terms are taken to be on 
boundary edges. 



4.1. Reliability: Testing with r e HK^^^. 
Lemma 7. Given K G Th, let 

{Ofork = 0, 
1 /2 
hK\Wh-Suh\\K + hi^ |||tr -kUhlWdK fork= 1, 
hKiWSahWK + hh-SuhWK) 
+h][\\\ltT *ahMdK + II [tr *UhMdK) for2<k<n. 

Let {a,u,p) be the weak solution to (3.6)-(3.8), let {(JhTUh,Ph) be the corresponding 
finite element solution having errors (eCT,e„,ep), and assume r S HA'^'^^{il) with 
llTllHAfc-i(o) < 1- Then 



|(e,,(T-n^T))-(d(T-n„T),e„)|< [ J2 ' 

\KeTh 



Proof. If fc = 0, then r is vacuous and so the term above disappears. We next 
consider the case 2 < k < n. Using Lemma 6, we write r = dip + z, where 
ip e H^K^-'^{n) and z € H^A^-'^in). Since UhT = dU^ip + U^z and dod^O, we 
have (i(T — n^r) = d{z — IIhz). Thus using the first line of (2.3) and the integration 



(4.3) 
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by parts formula (3.5) on each clement K G Th, we have 

-(eo-,r-n^r) + {d{T -n/ir),e„) = {(7h,T - H/it) - {d{T -XlhT),Uh) 

= X! ^'^''' '^('^ ^ -'^'''^) + z- T\hz)K - {d(z - nhz),Uh)K 
Ken 

Ken -^^^ 

+ {(Th - 5uh, z - Tlhz)K I tr (z - II/iz) A tr * u^. 

JdK 

Note next that tr (z — tt^z) is single- valued on an edge e = Kii^K2 , since z £ H^K^ 
and H-hZ G HK^. tr * u^ on the other hand is different depending on whether it 
is computed as a limit from Ki or from K2, and we use Jtr * u^l to denote its 
jump (|tr * u;i] = tr • u/i on 90). A similar observation holds for tr (cp — lih'^) 
and tr •cr/i. Let £h denote the set of faces (n — 1-dimensional subsimplices) in Th, 
and let -kQK denote the Hodge star on K^{dK) (with j determined by context). We 
then have using (3.3) and the fact that the Hodge star is an L2-isometry that 

y^ / tr ((^-H,i(p) A tr *crft = ^ (•ok tr(</?-H/i(^), |tr •cr,i]) 
(4.4) ^er. -^^^ eesu 

S Yl l|tr(v-n,,^)||aK||Itr •crJIIaK- 
KeTh 

Similarly manipulating the other boundary terms in (4.3) and employing (3.15) 
yields 

(cr/l,T - HhT) - (d{T - IlhT),Uh) 

< J2 V-i{K)[hK\\\z-nhz\\K + y-Ilh(p\\K) 
Ken 

+ /i^'/'dl tr (z - nhz)\\aK + II tr {(f - Ilh^)\\9K) 

^^•^) <( E V-i(Kry^'^{T. [/.^^(||^-H,^||^ + ||z-H,z|||,) 

Ken Ken 

+ /^K (II tr iv - Uh^WoK + II tr (z - H^z)i^)] ) 

< ( E v^Kf 

Ken 

Thus the proof is completed for the case 2 < k < n. 

For the case fc = 1 we have by definition that z = t E HA'^(n) = H^(il). Thus 
the proof proceeds as above but with terms involving 5ah and tr • (Jh omitted. D 

4.2. Reliability: Testing with v E HA''. In our next lemma we bound the term 
(/ — dah — Ph, V — H/jw) — {duh,d{v — H^^w)) from (2.26). Before doing so we note 
that ii" and ijJJ are always trivial, so in this case p = ph — 0- We however leave 
the harmonic term ph in our indicators even when k ~ n for the sake of consistency 
with the other cases. 



X 1/2 
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Lemma 8. Let K e Th, and assume that f G H A (K) for each K ^Th- Let 



(4.6) 77o(i^) = <^ 



hxWf -Vh - SduhWx + /i]/^|||tr -k duhlWoK fork = 0, 
hxiWf ~dah -Ph -Sduh\\K + \\S{f - dcr^ - ph)\\K) 
+h][\\\ltT *duh}\\9K + II [tr * if - dan - Ph)j\\dK) 
for 1 < fc < n — 1 
. \\.f-d(7h-ph\\Kfork = n. 



Under the above assumptions on the regularity of f and with all other definitions 
as in Lemma 7 above, we have for any v G HA''{n) with \\v\\HA''(n) ^ 1 



(4.7) {f-d<Jh-ph, 



- Hhv) - {duh, d{v - lihv)) < I Y, Vo{Kf ) 

\KeTh / 



1/2 



Proof. For k = n, the term (dw/i, d{v — Ilfiv)) is vacuous, and Galerkin orthogonality 

1 /9 

imphes that (/ - dah - Ph, v - Uhv) = (/ - d^h - Ph, v) < (Y^KeU "^oiK)^) 

For < fc < n — 1, noting that d(v — Ilfiv) = d{z — IlhZ+d{(p — Ilh(p)) = d{z — Ilhz) 
and integrating by parts yields 

{duh, d{v - Uhv)) = {duh,d{z - Uhz)) 

^ ' ' = \ {Sdufi, z — Hfiz) + tr (z — II/jz) A tr • dui^. 

Ken "^^^ 

For fc = both ip and ah are vacuous, so we may complete the proof by employing 
(4.8) and proceeding as in (4.4) and (4.5) to obtain 

{f-d(Jh -Ph,v- n^u) - {duh, d{v - Tlhv)) 

= Y^ {f - Ph - Sdufi, V - Uhv) tr (w - Uhv) A tr • duh 

Ken -^^^ 

< ( J2 ^oiK)' 

Ken Ken 

For l<fc<n — 1, we write v = dip + z as in Lemma 6 and employ (4.8) to find 

(/ - dah ~Ph,v - Xlhv) - {duh, d{v - \lhv)) 

= I (/ - dah - Ph, d{ip - Uhf))K 



(4.9) 



X 1/2 / ..-^ \ 1/2 



(4.10) 



+ y2 ^-f ^ '^'^^ ^Ph^ ^duh, z - Tlhz)K I tr (z - Whz) A tr • duh 

^Ken -^^^ 

^ [/] + [//]. 



The term // above may be manipulated as in (4.9) above in order to obtain 

(4.11) IIS[Y. h'K\\f-dah-ph-&duh\\']i + hK\\lirkduhl\\aK^^^ 

Ken 
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We now turn our attention to the term /. Integrating by parts while proceeding 
as in (4.4) and (4.5) yields 

I^ y^ (d{f - dah-Ph),^ -'nh(p)K 
Ken 



(4.12) + / tr (93 - Rhip) A tr * (/ - dcTh - Ph)- 

JdK 

1 /9 

<( E hl\\5{f-dau^PH)\?K + hK\\liT^{f-dan-pK)l\\lK) ■ 
Ken 

Combining (4.11) and (4.12) yields (4.7) for 1 < fc < n— 1, completing the proof. D 

We finally remark on an important feature of our estimators. The term /i^^ ||(5(/ — 
dah — Ph)\\K + hj^ llJtr * (/ — dah — Ph)}\\dK is in a sense undesirable because it 
requires higher regularity of / than merely / G i2- In particular, evaluation of the 
first term requires / e H*A'^{K) for each K e Th, and evaluation of the trace term 
requires tr * / G L2A{dK) for each K. (Note however that / is not included in 
the jump terms if / € H* A'' (il) .) Both relationships are implied by / G H^A'^{K), 
K E Th, so we simply make this assumption. 

In order to understand why such terms appear, note that the Hodge decompo- 
sition of / reads f = da + p + 5du. The first two terms da + p are approximated 
directly in L2 in the mixed method by da^ + Ph, while the latter term Sdu is only 
approximated weakly in a negative order norm (roughly speaking, in the space dual 
to HA*jS) in the mixed method. In our indicators, \\{da +p) — {dah +Ph)\\K is thus 
a naturally scaled and efficient residual for the mixed method, but \\Sdu — SduhHx 
is one Sobolev index too strong. The latter term should instead be multiplied by 
a factor of Iik in order to mimic a norm with Sobolev index —1, as in the term 
hK\\f - dah -Ph- Sduh\\K appearing in r]o. 

This "Hodge imbalance" implies that it is necessary to carry out a Hodge de- 
composition of / in order to obtain error indicators that are correctly scaled for 
all variables. When this decomposition is unavailable a priori, the Hodge decom- 
position must be carried out weakly in order to obtain a computable and reliable 
estimator in which the appropriate parts of the Hodge decomposition of / are scaled 
correctly. This is accomplished above. Since S{Sdu) — 0, hK\\S{f — dah —ph)\\K = 
hxiWSdCa + ScpWx). This scales roughly as a Sobolev norm with order —1 of 6dea 
and dcp, which in turn scales as the terms ||(ieo-|| + ||ep|j appearing in the original 
error we seek to bound. For an element face e G dfl, (3.9) along with tr -kp — on 
dfl imply that |tr • (/ — dah — Ph)l = tr • {da — dah — Ph) on dfl. Similarly, for 
an interior face e we have |tr * /] = |tr • {da — dah — Ph)l- 

If a partial Hodge decomposition of / is known, it is possible to redefine 770 

so that only f G L2 is required. If / = da + ip with ip = p + Sdu known a 

1/2 
priori, we may replace hK\\S{f - dah -Ph)\\K + hi^ ||[[tr -k{f -dah -Ph)l\\dK with 

1 /2 

hxWSphWK + \\da - dah\\K + hj^ \\lii -k ph\\\aK- li f = Q + Sdu with Q = da + p 
known a priori, we may instead replace this term with \\Q — dah—Ph\\K- We do not 
assume such a decomposition is known, since it if were one would likely decompose 
the Hodge Laplace problem into *8 and *B* problems, as described in [4]. 

We finally note that a similar situation occurs in residual-type estimates for the 
time-harmonic Maxwell problem curl curl u — w^u = /, where the elementwise indi- 
cators include a term /i/f|| div(/+a;^M/i)||x+/i^ ||(/+w^M/i)-n||aK (cf. [27]). There 
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however the assumption div f (z L2 is natural as it represents the charge density, 
and to our knowledge the connection of this term with the Hodge decomposition 
has not been previously explained. 

4.3. Reliability: Harmonic terms. Next we turn to bounding the terms in 
(2.26) related to harmonic forms. 

Lemma 9. Given qh G V^ and K G Th, let 

(4.13) V^{K,qh) ^ hK\\Sqh\\K + h]{^\\ltT ^ qhlWaK- 



Then if 1 < k < n and (f> E HA'' ^(O) with \\4>\\HA''-'-(n) ~ 1; 






(4.14) (<z,„ci(0-n^0))< ( J2 Vsj{K,qhf 



1/2 
r]^[K,qhyj . 
KeTh 

1/2 



Given an orthonormal basis {qi, ...,qM} for S)'^^, let fii — {TliKeT %(-^'9i)^) 
Then we additionally have 

1/2 
(4.15) ■ ■ '" ^\ 



g^v{^\^l)<^,:^[Y,^?^ 



Finally, if u^ G X , 

1 /2 

(4.16) ll^'84ll<( E ^^(^'"'^)') ■ 

Ken 

Proof. Let ip G iJ^A'^^^(fi) boundedly solve dip — d(p, as in (3.17) and preceding of 
Lemma 6. Employing (3.17), integrating by parts as in (4.8), and proceeding as in 
(4.9) immediately yields 



{qh,d{(j) - n,j0)) = E i^lhi 'P - ^iMk + tr ((y9 - Uhip) Atr-kqh 

T^^T- JdK 

(4-17) ,1/2 

Ken 

(4.15) immediately follows from (2.19) and (4.14), while (4.16) follows from (2.23) 
and (4.14). D 

4.4. Summary of reliability results. We summarize our reliability results in the 
following theorem. 

Theorem 10. Assume that fl C M" is a bounded Lipschitz domain of arbitrary 
topological type. Let < k < n. Let rj-i be as defined in Lemma 7, let tjq be as 
defined in Lemma 8, and let rjs) be as defined in Lemma 9. Let also {qi, ...qu} be 
an orthonormal basis for S)\^ and let /i be as in (4.15). Then 

l|ea||_H-Afc-i(o) + ||e„j|fl-Afc(o) + ||ep| 

^^•^^^ < ( E v-l{K)' + voiKr + m{PHfy^\^^\\u,\\. 

Ken 
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Let also uj^ be the projection of Uh onto 3^^- Then the term lJ,\\uh\\ in (4.18) may 
be replaced by 

(4.19) /x( E VfoiK,uif) +li'\\uh\\. 

Ken 

Proof. The four terms on the right hand side of (2.26) may be bounded by employ- 
ing Lemma 7, Lemma 8, Lemma 9, and once again Lemma 9, respectively. D 

5. Efficiency of a posteriori error estimators 

We consider efficiency of the various error indicators employed in §4 in turn. 
Before doing so, we provide some context for our proofs along with some basic 
technical tools. 

Efficiency results for residual-type a posteriori error estimators such as those 
we employ here are typically proved by using the "bubble function" technique 
of Verfiirth [25] (cf. [11, 6] for applications of this technique in mixed methods 
for the scalar Laplacian and electromagnetism) . Given K € 7h, ^ei bx be the 
bubble function of polynomial degree n + 1 obtained by multiplying the barycentric 
coordinates of K together and scaling so that maxa;^;^ 6/f = 1- Extending by 
outside of K yields bx <E M^^(rj) with supp(5i<-) = K. Similarly, given an n — 1- 
dimensional face e — Ki D K2, where Ki, K2 G Th and K2 is void if e C dQ, we 
obtain an edge bubble function be defined on Ki , K2 by multiplying together the 
corresponding barycentric coordinates except that corresponding to e and scaling 
so that Taaxjii &e = 1- 

Given a polynomial form v of arbitrary but uniformly bounded degree defined 
on cither K G Th or a. face e C K G Th, 

(5.1) \\v\\k ^\\\/lrKv\\K, ||u||e ^ llV^^lle- 

Also, given a polynomial k-ioira v defined on a face e — Ki U ^2, we wish to define 
a polynomial extension Xv of v to Ki U K2. First extend v in the natural fashion 
to the plane containing e. We then extend v io Ki, i = 1,2 by taking Xv to be 
constant in the direction normal to e. Shape regularity implies that e, Ki, and K2 
are all contained in a ball having diameter equivalent to hx ~ hxi — hx^ , so that 
an elementary computation involving inverse inequalities yields 

(5-2) \\Xv\\l^(K^UK2) S hK^ML2(e)- 

5.1. Efficiency of 77_i. We first consider the error indicator ?7_i. 
Lemma 11. Let K E Th- Then for 1 < k < n, 

(5.3) ?7-i(-ft^) < ||eo-|U^ + ||e„||^^. 

Proof. We begin with the term hK\\<yh — ^UhWn- Let ip = bKi^^h — 5uh) <E HA'^^^; 
note that It dKi' — 0. Employing (5.1), the first line of (2.3), and the integration- 
by-parts formula (3.5), we obtain 

2 



0-/1 - 5uh\\K - (o-ft - Suh,ip) = {(Th - cr, -0) + {dip, u - Uh) 
\\kM\k + \\eu\\K\\dip\\K- 



^^■^^ <ll^-ll-||^lk + ||e„||K||rf^| 



A POSTERIORI FEEC 23 

Employing an inverse inequality and 6_ff < 1 yields 11^-011^^ < h^ \\'4'\\k ^ h^ llo'ft ^ 
dufiWii- Multiplying (5.4) through by /i/f/||(T;i — SufiW^ while noting that hj^ < 1 
yields 

(5.5) hK\\a-h - Suh\\K < ||ecr|k + ||e«||K- 

Let now fc > 2. Recall that Sa ~ SSu = 0. Thus with ijj = hx^Uh, we have 

||(5(T/j||^ ~ ((5crft,V) = ((5(crft -cr),V') = (c^h -cr, rfV') 

(5.6) _ _ 

< hj^ \\ea\\K\\i>\\K < hj^ \\ea\\K\\5(yh\\K- 

Multiplying through by Hk /\\5(Th\\K yields 

(5.7) hK\\5<Th\\K < \\ea\\K- 

We now consider edge terms. Note from (3.7) and the fact that u € H*K^{^) 
(since 5u = a € L2h^~'^{^)) that we always have |tr *u] = (suitably interpreted 
in i7-i/2). Let £h^ e = Kif] K2, where ii'2 = if e C 90. Jtr •u^] e A'=(e), so 
we let -0 e A"~-'^~'^(e) satisfy ^^ — |tr *Wft,]. The definition of* implies that ■(/; is a 
polynomial form because |tr ^w^l is. Note also that multiplication by h^. commutes 
with tr and 7k- since both are linear operations, so that h^i^tp ~ •kibeip) — *tr (feeXi/j). 
Employing the polynomial extension Xip defined in (5.2) and surrounding along with 
the second relationship in (3.3) thus yields 

|||tr *Uhj\\l - (foe*V', Itr -kUhDe 

(5.8) f 

= (*tr(6eXi/')Jtr *Uh])e = / tr (beXi/-) A [[tr • m/i]. 

J e 

Employing the integration-by-parts formula (3.5) individually on Ki and K2 yields 

(5.9) / tr(6eX^) A [tr -kUhj == {d{beX-4,),Uh)KiuK2 - {beX^,^hUh)KiuK2- 

J e 

Here Sh is 6 computed elementwisc, which is necessary because Uh ^ H* A^ glob- 
ally. Also, heXii ^ -ffA'^^^(O) with support in Ki U K2. Inserting the rela- 
tionship {(Jjb^Xip) ~ {d{he,Xtp)i'^) ~ into (5.9) and using an inverse inequal- 
ity |lrf(&eXi/')ll-ff ^ ^li W^sX^Wk^ (5.2), and the Hodge star isometry relationship 
\\%l>\\e ~ |||tr •Uftllle then yields 

jlltr *M/il||e ^ (cr, beXi/-) - {beX->p,^hUh) + {d{beXip)iUh -u) 

= (ecr, 6eX^) - {beXi(:,<^h - ShUh) + (d{beXi:)i e„) 

^ ll&eXl/'ll^lUKgdlecrllAriUKg + ^ j^ 1 1 e„ 1 1 _ff ^y Ka + \\<^h — (^hUh\\ K1UK2) 

< ft.]/^|||tr ■kUhl\\e{\\ea\\KiUK2 +^K^||e„||KiUK2 + \Wh - Suh\\KiUK2)- 

Multiplying both (5.10) through by h^^ /llltr *u/i]||e and employing (5.5) thus 
finally yields 

(5.11) h]/^\\ltT ^UhlWe < \\eu\\K,uK2 + \\eAK,uK2- 
A similar computation yields 

(5.12) h][^\\lti^a,Me<\\ejK,uK2, 

thus completing the proof of Lemma 11. D 
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5.2. Efficiency of rjo. We next consider the various error indicators 770 in Lemma 
8. First we define three types of data oscillation. First, 

(5.13) 08c{K)^hK\\.f-P.f\\K, 

where Pf is the L2 projection of / onto a space of polynomial k-iorms of fixed but 
arbitrary degree. Note that Pf is in general globally discontinuous. We do not 
specify the space further, since it is only necessary that it be finite dimensional in 
order to allow the use of inverse inequalities. Similarly, we define the edge oscillation 

(5.14) osc(ai^) - h]l^\\ltT * (/ - PfML^OK)- 
Finally, we define 

(5.15) OSCsiK) - hKim - Pf)\\L,(K)- 

For a mesh subdomain lj of ft, let osc((x;) = (X^tcw osc(iir)^) and similarly 
for osc s . The last two oscillation notions measure oscillation of da only, since the 
Hodge decomposition yields Jtr * /] = |tr Tk-dcr] and Sf — 5da. 

Lemma 12. Let K G Th, and consider the error indicators defined in Lemma 8. 
For k — 0, we have 

(5.16) T]o{K) <\\eu\\H,ujK + osc (ujk) 
When k — n, 

(5.17) m{K)<\\de„\\K + \\ep\\K- 
Forl<k<n-l, 

Vo{K) <||e„||H,^^, + lle^llH,^^ + |iep|l^^ 

+ osc{ljk) + oscs{luk) + osc{dK). 

Proof. For the case k = n, (5.17) follows trivially from the Hodge decomposition 
f — da + p and the triangle inequality. 

For the case < fc < n — 1, let V' = bx{Pf — dah — Ph ^ Sdui^). Then 

\\P.f -dah -Ph -5duh\\]i ~ {Pf -dah - Ph - 5duh,ip)K 

= {Pf - f, V'>K + (/ - dah -Ph- 5duh,-tp)K- 

Employing the Hodge decomposition / = da + p + Sdu and then integrating by 
parts while recalling that bx and thus ip vanishes on dK yields 

(./ - dah -ph- Sduh, -iP)k = {dca- + e^ + Sdeu, ^)k 

^ {e„,Sip)K + {ep,rp)K + {deu,d'ip)K- 

Collecting (5.19) and (5.20) and then employing the inverse inequality ||(ii/;||if + 
ll'^V'll-ftr ^ h~^ II "01 1 AT I multiplying the result through by Hk, and dividing through 
by \\Pf-dah-ph-Sduh\\K after recalling that WiPWk < \\P f - dah - Ph - Sduh\\K 
yields 

(5.21) hxllPf -dah - Ph - Sduh\\K < \\eAK + \\deu\\K + hK\\ep\\K + osciK). 

Employing the triangle inequality completes the proof that /ik||/ — dah — Ph — 
Sduh\\K is bounded by the right hand side of (5.18) when 1 < fc < n — 1, or by 
the right hand side of (5.16) when k = after noting that in this case p = ph is a 
constant and recalling that cr — ct^ is vacuous. 
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We next consider the term /i^ |||tr ^(iw/j]||^x in (4.6). Manipulations similar 
to those in the previous subsection yield 

1 /2 

(5.22) hj^ \\ltr-kduh}\\e<\\deu\\KiuK2+hK\\Sdeu\\KiuK2- 

Employing the Hodge decomposition / = da + p + 6du yields Sde^ — (f — dat — 
Ph — Sdufi) — dca — Cp. Thus 

1 /2 

f5 23l ^^ llltr *rfuj||e < ||de„||xiuA:2 

+ hKiWf - dah - Ph - Sduh\\KiuK2 + \\dea\\KiuK2 + \\ep\\KiuK2)- 

Employing (5.21) on Ki and K2 individually completes the proof that /i^ ||[[tr * 
dw/iJIle is bounded by the right hand side of (5.18) in the case 1 < fc < n -- 1 and 
of (5.16) when fc = 0. 

We now consider the term hK\\S{f — da^ — Ph)\\K- First note that hK\\5{f — 

d<Jh-Ph)\\K < oscs(K) + hK\\S{Pf-d(Th-Ph)\\K- Letting?/; = bK5{Pf -dcTh-Ph) 
and recalling the identities 5f — Sda and dp ~ 0, we integrate by parts to compute 

\\S{Pf - dah - Ph)\\l c^ {SiPf - dah - Ph),i^) 

= {5{Pf-f),yj) + {5{de, + ep),i:) 

< h]^^oscs{K)\\tp\\K + \{de„ + ep,d^)| 

< h]^^{oscs{K) + \\de„\\K + ||ep||x)||'0lk- 

Further elementary manipulations as in (5.19) and following complete the proof 
that hK\\S{f — d(Jh — Ph)\\K is bounded by the right hand side of (5.18). 

1 /2 

We finally turn to the edge term hj^ |||tr * (/ — duu — Ph)\\\aK- Note first that 
|tr •(p+Jdu)] = on all element faces e. On interior faces this is a result of the fact 
that p + 5du E H* K^ , while for boundary edges this is a result of (3.9) along with 
the definition oi ^^ . Thus |tr •/] = |tr -kda}. Setting TtrV = [tr ^{Pf-duh-Ph)} 
and letting x^ be the polynomial extension of ij: as above, we compute for a face 
e^ Kif]K2 that 

II Itr * (P/ - dGh - Ph)\ \\l ^ {be * ^, Itr * (Pf - dah - Ph)l) 
< h^'/\sc{dK)\\i;\\e + \{be^, Itr * (de, -pj)| 



(5.25) 

+ {beX-4>,S{dea- -ph))KiUK2\- 



1 /o 

= /l^ OSc(aii:)||V'||e + \{d{beX4,),de^ -Ph)KiUK2 



Next note that {d{beXip),p) — 0, so that {d{bex4') , ^Ph) = {d(beXip) , Sp) ■ Using an 
inverse inequality and (5.2) then yields 

llltr i^iPf -dah-ph)j\\l < hj^'"^ osc{dK) + \\ep\\K^uK2 + \\dea\\KiuK2 
(5.26) ^ 

+ hK\\S{f -d<Th-ph)\\KiUK2 Me- 

1/2 
Further elementary manipulations as above complete the proof that /i^ || Jtr • (/ — 

d<Jh — Ph)l\\e is bounded by the right hand side of (5.18). 

D 
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5.3. Efficiency of harmonic indicators. We finally state efficiency results for 
the various harmonic terms. 

In this section we prove efficiency of the individual harmonic terms appearing 
in Lemma 9. As we discuss more thoroughly below, however, we do not obtain 
efficiency of all of the terms that we originally sought to bound. 

Lemma 13. Let Vh <E Vj^ . Then 

(5.27) m{K.Vh)<\\P^Vh\\u,^- 
In particular, we have for uj^, qi G i^^, and ph 

(5.28) m{K,ui) < \\P<3uUU, 

(5.29) V^{K,q^) < \\P<sq^L^^\\q^-PmL^, 

(5.30) V^iK,Ph) < \\ep\U. 

Thus with II and fii as in Lemma 9, 

(5.31) iJi<gM^\^h) 

Proof. The proof of (5.27) is a straightforward application of the bubble function 
techniques used in the previous subsections. (5.28) and (5.29) are special cases of 
(5.27), while (5.30) may be proved similarly. Finally, summing (5.29) in H.2 over 
K (z Th while employing the finite overlap of the patches lok (which is a standard 
consequence of shape regularity) implies that 

(5.32) l^z^hi- PsjQiWiJKi 

which yields (5.31) when summed over 1 < i < M. D 

Remark 1. Lemma 13 gives efficiency results for the terms in our a posteriori 
bounds for gap(i5'^, Jojj) and for HPigM^H, but not for the quantity ||Pf,M;i|| itself 
that we originally sought to bound. More generally, we have not bounded all of 
the harmonic terms (4.18) and (4.19) by the error on the left hand side of (4.18) as 
would be ideal. The offending terms arc due to the nonconforming nature of our 
method which arises from the fact that Sj'j^ ^ S)^ . Establishment of efficiency of 
reliable estimators for this harmonic nonconformity error remains an open problem. 

6. Examples 

In this section we translate our results into standard notation for a posteriori 
error estimators in the context of the canonical three-dimensional Hodge-de Rham 
Laplace operators. Below we always assume that n = 3. 

6.1. The Neumann Laplacian. When /c = 0, cr and the first equation in (2.3)) 
are vacuous. Also, V^-^ = V-^ =%,V^ = V^ = H^{n), d = V, and (5 = - div. In 
addition, p = Ph ~ /q /j ^^'^ ^'^^'^ — ^Au. The weak mixed problem (2.3) reduces 
to the standard weak form of the Laplacian and naturally enforces homogeneous 
Neumann boundary conditions. In Lemma 7 and Lemma 9 we have ry-i = and 
/^ = ll^s^ft II — Oj respectively. Thus r/i is thus the only nontrivial indicator for 
this problem, and it reduces to the standard indicator ri{K) from (1.3). Thus we 
recover standard results for the Neumann Laplacian. 



A POSTERIORI FEEC 27 

6.2. The vector Laplacian: fc = 1. For fc = 1 and fc = 2, the Hodge Lapla- 
cian corresponds to the vector Laplace operator curl curl —V div, but with different 
boundary conditions. For A: = 1, m G iJ(curl), u = — divM € H^ , and the boundary 
conditions are u ■ n — 0, curlw x n = on dQ,. S^^ consists of vector functions p 
satisfying curlp ~ 0, divp = in f7 and p • n = on 951. From (4.1), 

(6.1) r;_i(ii') = ft,/f||o-/j + divM/j||K + /i^ ||[[uh • n]||ax- 

Here n is a unit normal on dK. From (4.6) we find 

rin(K) = Zixdl/- '^(^h -Ph -curl curl M/i II K + || div(/ - Vdh -ph)\\K) 

+ h]l\mcnTluh)tMdK + \m-Vah-ph)-nl\\9K), 
where a subscript t denotes a tangential component. Finally, in Lemma 9 we have 

(6.3) m{K,qh) = hxWdWqhWK + hj^ II kh • "lllaif- 

6.3. The vector Laplacian: fc = 2. In the case k — 2 the mixed form of the 
vector Laplacian yields a ~ curlw, u G _ff(div), and u x n — 0, divu = on 90. In 
addition, Jo^ consists of vector functions p satisfying curlp = 0, divp = in O and 
p X n ^ Pt ~ on dil. We then have from (4.1) that 

V-i{K) = hxiW divcr^llx + ||cr/i - curlu^H/f) 

+ h]['m<J,,-nj\\oK + \\lu^,,tl\\aK)- 
From (4.6) wc have 

rio{K) = hKiWf - cuTl(Th -Ph + Vd[vuh\\K + \\ curl(/ - curlcr/i -ph)\\K) 

h]l\\\ldivuhMdK + mf-cnTlah-pH)tMdK). 
Finally, in Lemma 9 we have 

1 /2 

(6.6) r]^{K,qh) ^ hK\\cuT\qh\\K + hj^ II fo, til lax- 

6.4. Mixed form of the Dirichlet Laplacian. For fc = 3, (2.3) is a standard 
mixed method for the Dirichlet Laplacian — Au = in il, m = on dfl, and 
a = -Vm. d^ = div, d^ is vacuous, S)^ = Sj^ = 0, V''^'^ = iJ(div), and V'' = i2- 
Taking ah and Uh now to be proxy vector fields for a^ and Uh, we have in (4.1) 
that 5(7h = curler^, 5uh = — Vm/i, tr •cr^ = a^.t (i-c, the tangential component of 
ah), and ti *Uh = Uh- 

(6.7) f]-i{K) - hKiW cuTlahU + hh + "^UhU) + 4^'(IIK Jlk + IIKlIk)- 
In addition, (4.6) yields 

(6.8) r;o(ir) = ||/-diva^||;f. 

The "harmonic estimators" in Lemma 9 are all vacuous in this case. Combining 
Theorem 10 with the corresponding efficiency bounds of §5 thus yields 

(6.9) ||en||L.(^2) + l|e.||H(div;^2) ^ ( ^ ^-i W + m(K)^f^'- 

Ken 
In contrast to the vector Laplacian, many authors have proved a posteriori error 
estimates for the mixed form of Poisson's problem, so we compare our results with 
existing ones. We focus mainly on two early works bounding a posteriori the natural 
mixed variational norm iJ(div) x L2- In [8], Braess and Verfiirth prove a posteriori 
estimates for ||eCT||ff(div;n) + l|e«ll) as we do here, but their estimates are only valid 
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under a saturation assumption (which is not a posteriori verifiable) and are not 

efficient. Sahent to our discussion is their observation on pp. 2440-2441 that 

the traces of iJ(div) test functions he only in H^^'"^. This prevented them from 

employing the mixed variational form in a straightforward way, that is, using an 

inf-sup condition in order to test with functions in iJ(div) x L2. Doing so using 

their techniques would have led to a duality relationship between traces lying in 

incompatible spaces, or more particularly, between traces lying in H~^'^ and some 

space less regular than H^/^. Following ideas used in [11] in the context of the 

mixed scalar Laplacian and developed more fully in [24] for Maxwell's equations, 

we insert the essential additional step of first taking the Hodge decomposition of 

test functions. Only the regular [H^) portion of the test function is then integrated 

by parts, thus avoiding trace regularity issues. Note finally that the elementwise 

1/2 

indicators of [8] are of the form || d\\(Th—f\\K+\Wh+'^Uh\\K+hj^ \\\uh\\\dK, which 

includes our indicator 779 and parts of our indicator ?7-i. However, the jump term 

^K \\\'^h\\\dK is scaled too strongly (by hj^ ' instead of hj^ in our estimator), 

and the resulting bounds are thus not efficient; cf. (4.20) of [8]. 

In [11] Carstensen provided a posteriori estimators for the natural iJ(div) x 

L2 norm which are equivalent to the actual error as in (6.9). In our notation, 

Carstensen's elementwise indicators have the form ||/ — diver ;j||j^; + /ixjj curl(j/j||if + 

1 /2 

ft-if min^j^gL^ \\ah + Vuh||x + hj^ \\l(Th,t\\\dK- Here Lu is an appropriate space of 

1 /2 

piecewise polynomials. Thus our terms /iK||crh+VM/i||i<-+/i^ || [[uh]||aK are replaced 
in Carstensen's work by Hk min^^gi^ \\ah + Vw^||i^:, and our estimators are other- 
wise the same. However, Carstensen's results were proved only under the restrictive 
assumption that 51 is convex, which we avoid. [11] also makes use of a Helmholtz 
(Hodge) decomposition, but a commuting quasi-interpolant was not available at the 
time and thus full usage of the Hodge decomposition was not possible. 

Most works on a posteriori error estimation for mixed methods subsequent to 
[11] have focused on measuring the error in other norms, e.g., ||ecr||L2 (cf- [19, 26]). 
One essential reason for this is that the i7(div) x L2 norm includes the term ||/ — 
divcr^ll which directly approximates the data / and which can thus be trivially 
computed a posteriori. Thus while the _ff(div) x L2 norm is natural to consider 
from the standpoint of the mixed variational formulation, it is perhaps not the most 
important error measure in practical settings. Even with this caveat, our estimators 
for mixed methods for the Dirichlet Laplacian seem to be the first estimators that 
are directly proved to be reliable and efficient for the natural mixed variational 
norm under reasonably broad assumptions on the domain geometry. 
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